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Abstract 



> 

A model many-body Hamiltonian describing an heterogenous system of 

o 



paired protons and paired neutrons and interacting among themselves through 
monopole particle-hole and monopole particle-particle interactions is used to 
study the double beta decay of Fermi type. The states are described by time 
dependent approaches choosing as trial functions coherent states of the sym- 
metry groups underlying the model Hamiltonian. One formalism, VP1, is fully 

X 

$— 1 ' equivalent with the standard pnQRPA and therefore fails at a critical value of 

the particle-particle interaction strength while another one, VP2, corresponds 
to a two step BCS treatment, i.e. the proton quasiparticles are paired with 
the neutron quasiparticles. In this way a harmonic description for the double 
beta transition amplitude is provided for any strength of the particle-particle 
interaction. The approximation quality is judged by comparing the actual 
results with the exact result as well as with those corresponding to various 
truncations of the boson expanded Hamiltonian and transition operator. Fi- 
nally it is shown that the dynamic ground states provided by VP1 and VP2 
are reasonable well approximated by solutions of a variational principle. This 



remark constitutes a step forward finding an approach where the RPA ground 
state is a solution of a variational principle equation. 



I. INTRODUCTION 

One of the most exciting subject of the modern nuclear physics is the double beta decay. 
The reason is that from this field one expects an answer to the question of whether neutrino 
is a Majorana or a Dirac particle. The process may take place through one of two channels, 
Ovfl/3 and 2vf3j3. The answer to the above mentioned question might come from the discov- 
ery of the first process which if it exists requires that the neutrino is a Majorana particle 
and which should be a rare one |l]-f§. In order to make predictions of the neutrino mass and 
the right handedness of its electroweak interactions one should use reliable nuclear matrix 
elements. However there is no stringent test for these matrix elements. Fortunately for the 
double beta decay with two anti-neutrinos in the final state 2v/3@, one uses similar nuclear 
matrix elements and moreover for this process plenty of data exist. Therefore the idea of 
using for the neutrino-less double beta decay the same many-body approach and the same 
NN force which describe realistically the two neutrino double beta decay, was adopted by 
most of the groups involved in such studies. Actually this is the reason why so many theo- 
reticians focused their attention in explaining the features of the experimentally measured 
2u(3(3 process |]7|-[T3|. The formalism which predicts decay rates closest to the experimen- 



tal data is the proton-neutron quasiparticle random phase approximation (pnQRPA) which 
includes the particle-particle (pp) two body interaction. Without this interaction the pre- 
dicted transition amplitude is too large as compared to the existent data. Including the pp 
interaction and considering its strength, g pp , as a free parameter one obtains the following 
behavior for the transition amplitude. For a large interval of g pp , starting from zero, there 
exists a plateau followed by a quick decrease. The amplitude vanishes at a certain g pp and 
shortly after this point one reaches, with increasing g pp , the critical value where the pnQRPA 
breaks down. On the decreasing part of the 2vj3(3 transition amplitude as a function of g pp 
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one finds agreement with experiment. The drawback of this description is that in this area 
of g pp , the pnQRPA is not a good description and therefore the results are not stable to 
adding higher RPA correlations. Many attempts have been made to stabilize the ground 
state in this region of the particle-particle interaction. The new methods tried either to keep 
the pnQRPA boson picture but include higher correlations through boson expansion JTjj-jTjl 
techniques or via self-consistent procedures |20fl , or to re-normalize the pnQRPA phonon 



|2TH24|. These methods improved the description of the data j25| although they have also 
drawbacks and lack sometimes of consistency |24],^,[^| . 



In a previous paper |26| we advanced the idea that the ground state may be stabilized by 
introducing the particle-particle interaction first in the mean field and then to the pnQRPA 
process. A way to do that is even possible at the level of independent quasiparticle repre- 
sentation. Indeed, we noticed in a schematic single j model |29| that one term of the two 



body quasiparticle interaction is a quasiparticle proton-neutron pairing interaction and its 
strength is negative. Therefore one could define a Cooper pair out of one proton quasiparticle 
and one neutron quasiparticle. This was worked out through a time dependent variational 
principle. We chose alternatively four distinct trial functions. Two of these define different 
RPA approaches. One is identical to the standard pnQRPA and works for small strengths 
of the particle-particle interaction, namely before the critical point, whereas the other one 
provides a harmonic picture beyond the critical value of the interaction strength, i.e. in the 
region where the standard pnQRPA does not work at all. Moreover in this region a boson 
expansion procedure might be defined and through diagonalization the exact result can be 
recovered. 

In the present paper we continue the study started in the previous paper, by focusing on 
the following new features. Here we calculate the double beta decay transition amplitude 
by using the description in the interval of the particle-particle interaction forbidden for the 
pnQRPA approach. We hope that by matching the approaches in the two complementary 
intervals one could provide a unified description of the process for small as well for large 
strengths of the particle-particle interaction. When one uses the boson expansion for the 



model Hamiltonian one aims at testing the convergence properties by comparing the cor- 
responding results with the exact ones. One should mention that this is possible only for 
a single j model Hamiltonian involving a monopole two body interactions. Formulating it 
for the double beta decay Fermi transition the boson expansion of the transition operator 
emitted in ref. [TTJJ is a Schwinger type boson expansion since it uses two different bosons 
in order to satisfy the condition of mapping two algebras. This is a good test for the boson 
representation which was previously used. Another aspect which is treated here refers to the 
structure of the vacuum states provided by the standard pnQRPA and the approach which 
works beyond the critical particle-particle interaction strength. We want to check whether 
these two vacua might be described in a reasonable approximative fashion by static ground 
states yielded by a variational principle with a suitable trial function. 

This project is achieved according to the following plan. In Section 2, a brief review of 
the results obtained in the previous paper is given. These results are used here for a self- 
consistent presentation. In Section 3, we compare the results for double beta decay transition 
amplitude obtained in the two intervals of the particle-particle interaction, i.e. before and 
after its critical value, with different methods. In Section 4 the Holstein-Primakoff boson 
representation is used for the model Hamiltonian. By diagonalizing the full Hamiltonian, 
the first order approach and the results for the limit Q — > oo and by comparing the results 
with the exact ones and those obtained within the harmonic picture, one could judge the 
convergence quality of the expansion used. The Schwinger type boson representation is here 
also applied. The possible contribution of pairing vibrational states to the enhancement of 
the decay rate is discussed. The decay to the double phonon pairing vibration is explicitly 
analyzed. In Section 5 we study a possible relationship between the vacuum states of the 
standard pnQRPA, the second order BCS approach, and the static ground state provided by 
classical variational equations obtained with two different trial functions. Final conclusions 
are drawn in Section 6. 
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II. BRIEF REVIEW OF SOME PREVIOUS RESULTS 



Since the present paper is continuing the study we started in previous publications we 
give here a brief summary of the main ingredients used there |26[. In this way we fix 
the notation and the conventions. Also the general frame of our present results will be 
better emphasized. The object of our study is a system of protons and neutrons moving in a 
spherical shell model mean field, interacting among themselves through pairing, particle-hole 
and particle-particle monopole two body interaction. The associated Hamiltonian reads: 

H = € p C prn C prn + e n C nm C nm 
m m 

pip G n pfp 

2 v p 2 

+ 2 X /3-/3 + ~ 2*iP-P + , (2.1) 
where the following notations have been used: 
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m 
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p- 


= Tc l $ ■ 


P 


f = (p-y 



(2.2) 

The standard notations for the creation and annihilation operators are used. For the sake 
of simplicity we consider here the case of a single j shell. The extension to the multi-shell 
picture is straightforward. However for the purposes of the present study the multi-shell 
calculations are not necessary. The shell j is self-understood in our notations. To specify 
the isospin of the particle occupying the shell we use the indices p for proton and n for 
neutrons. The time reversed state is denoted by \ fm) = (— l) J_m |r, — m). Performing the 
Bogoliubov-Valatin transformation, the model Hamiltonian is transformed in a many body 
quasiparticle operator which after ignoring the quasiparticle scattering terms and, for the 
sake of simplicity, equating the proton and neutron single particle energies, has the form: 



H = e{N p + N n ) + XiAl n A pn + A^^-dJn + A pn A pn ), (2.3) 
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where N p and N n are the proton and neutron quasiparticle number operators, respectively, 
while A+ n , A pn stand for the proton-neutron pairing operators built up with the quasiparticle 
operators aj m ,4 m : 

4n = K a n]oo, A P n = (AlJ. (2.4) 

The coefficients e, Ai,A2 are functions of the parameters defining the Hamiltonian (2.1) as 
well as of the coefficients U, V determining the quasiparticle representation. The Hamilto- 
nian (2.3) is a quadratic expression in the generators of the SU(2) algebra 

S = ^[N p + N n -2Q], S + = V2QAl n , S- = V2QA pn , Q=?l±l. (2.5) 



and therefore it is exactly solvable. This feature has the great advantage that one could test 
the many body approximations adopted in various formalisms. This virtue made it very 
attractive for many theoreticians who used it widely to advance some new hypothesis for 
treating the proton- neutron two body interaction |29|-|3T|| . 

One may argue that restricting the model space to a single j level one misses the effect 
of spin-orbit interaction. However from our earlier studies |3"2],|3"3]] , we now that the spin 



orbit coupling is important for the Gamow- Teller double beta transition, where most of the 
strength is carried by spin-flip configurations, but not for the Fermi transition which we are 
studying in the present paper. 

If one assumes that the monopole two quasiparticle operators and their hermitian conju- 
gate satisfy bosonic like commutation relations than a harmonic excitation operator might 
be defined. Therefore the Weyl algebra, defined by the operators A', A, I (here I denotes 
the unity operator), can be used to describe the eigenstates of H. When the terms with 
four quasiparticles in H are dominant and the quasiboson approximation is adopted, for the 
above mentioned commutators, one could diagonalize first the anharmonic part of H and 
then treat the term expressing the quasiparticle total number operator. Up to an additive 
constant, the anharmonic part of H (four quasiparticle terms) is a linear combination of the 
generators of the group SU(1, 1), i.e. A +2 ,A 2 ,2 + AA + A. Therefore this symmetry might 
be useful for describing the eigenstates of H. 
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The model Hamiltonian (2.1) was treated, in ref. p6 |, [] within a dependent variational 
principle formalism: 





5 j^\H-i— \i>)dt' = Q. 



(2.6) 



by taking as trial functions the coherent states of the above mentioned symmetry groups. 
Since the coherent state for SU(2) and for the Weyl group are formally identical with the 
difference that in the first case the operators A' , A satisfy exact commutation relations, while 
in the later case they are considered as quasibosons, we complete our study with a fourth 
wave function which is similar to the SU(1,1) coherent state but the operators involved 
satisfy exact commutation relations. Concluding, the following trial functions have been 
alternatively considered for the variational state |26,27|: 



l^iCz, z*)) = exp{zA\ n - z*A pn }\0) q , \A pn , A\J = 1, (VP1) 



\ip2(z, «*)) = exp{zA\ n - z*A pn }\0) q , \A pn , A\ n 



\fo(z, z*)) = exp(^J n 2 - z*A pn 2 )\0) q , \A pn , At 



pn 



1 

~ 2fi 
1, (VP3) 



(N p + N n ), (VP2) 



\Mz,z*))=Af 4 e z < 2 \0) q , [ Apn ,Al n } = l-^{N p + N n ) ) (VP4) 



(2.7) 
(2.8) 
(2.9) 
(2.10) 



In each case, z is a complex function of time and z* denotes its complex conjugate. |0) 9 
is the quasiparticle vacuum state. The formalism corresponding to the trial function |^) 
with the commutation relations specified above, will be hereafter called VPk. For each of 
the four cases we solved the static equations of motion and then found the harmonic mode 
describing small oscillations around the static ground state. Several properties such as the 
behavior of energies with respect to the strength of the particle-particle interaction, as the 
ground state correlations, as the single beta transition amplitude, as the Ikeda sum rule and 
also various quantization procedures of the classical equations of motion have been analyzed. 



For a better presentation the final results will be considered, as in ref. p6| , p9| , as function 
of the re-scaled strengths of particle-hole and particle-particle interactions: 



1 Throughout this paper we adopt units where h=l 
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x' = 2n x , k' = 2n X i. (2.11) 

For example the energies of the first excited states in mother and daughter nuclei, predicted 
by VP1 and VP2, vary with k' as in Fig 1. while those given by VP3 and VP4 as shown in 
Fig 2. 

Here we complete this analysis with some new properties of the proton-neutron interact- 
ing system. 



III. DOUBLE BETA DECAY 

The double beta decay with two anti-neutrinos in the final state is considered to take 
place by two consecutive single (3~ decays. The intermediate state reached after the first 
beta decay, consists of an odd-odd nucleus in a pn excited state, one electron and one anti- 
neutrino. If in the intermediate state, the total lepton energy is approximated by the sum 
of the electron mass and half of the Q- value of the double beta decay process (AE), the 
inverse of the process half-life can be factorized as follows : 

(T^)" 1 = F | Mp | 2 , (3.1) 

where F is a lepton phase integral while the second factor is determined by the states 
characterizing the nuclei involved in the process and has the expression: 

M y ^ ^(0 + |/g + |0^)^(0^|0^) drf (0^|/3 t |0 + )d , , 

F U> E k + AE ' 1 ' ' 

where the transition operators is defined by eq (2.2). The initial and final states are the 
ground states of mother (|0^) m ) and daughter (lO^) nuclei. These states are vacuum states 
for the pnQRPA phonon operators built on the top of the corresponding BCS states which 
are the static ground states. The matrix elements describing the first and the second legs 
of the process have been calculated in a previous publication based on the semi-classical 
approaches (see eqs. 6.25 and 6.32 from ref. |26[| ) VP1 and VP2. In both cases there is 
only one excited states and therefore the summation in eq. (3.2) can be omitted. The 
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overlap matrix element involved in Eq. (3.2) is function of the forward (X) and backward 
(Y) amplitudes of the harmonic boson and has the expression: 

m(0 ] f |0+)<2 = X m X d — Y m Y d , (3.3) 

where the indices m and d specify that the corresponding states describe the mother and 
daughter nuclei, respectively. The energy denominator in Eq (3.2) is: 

£i = u m , (3.4) 

where u m , is the energy of the harmonic mode in the mother nucleus. The results for the 
double beta transition amplitude obtained within the VP1 and VP2 approaches are plotted 
in Fig. 3 a) as function of the proton- neutron pairing strength k'. In the VP1 approach the 
energy shift AE is taken equal to 1 MeV. The same AE is considered for the first part of 
the k' interval in Fig. 3 b). For the strength of the proton- neutron pairing k' > l.3MeV 
this energy shift is corrected by half the difference between the static ground state energies 
in the mother and daughter nuclei. From there we see that the transition amplitude for 
small (k' < 0.8MeV) and large (k' > 1.5MeV) k' is almost constant. Moreover there is a 
finite interval for k' where the amplitude is not defined. The justification for such a window 
can be seen in Fig 1. Indeed from the lower panel it is clear that the VP1 approach yields 
a transition amplitude which is not defined beyond k' = 1.05 since the harmonic mode in 
the daughter nucleus collapses. In the upper panel of Fig.l it is shown that the SU(2) mode 
in the mother nucleus has a non- vanishing energy only for values of k' larger than 1.2 for 
X ' = 0. and 1.3 for X ' = 0.5. 

IV. BOSON EXPANSION 

For the VP2 formalism several classical canonical coordinates have been found. To each 
set of coordinates corresponds a specific quantization scheme. Through the quantization 
procedure, the classical energy function is transformed to an operatorial function of boson 



9 



operators. In what follows we shall use two boson representations of the model Hamiltonian. 
In each case the corresponding eigenstates are used to evaluate, by means of eq. (3.2) , the 
transition amplitude for the double beta Fermi transition. 



A. Holstein Primakoff boson expansion 



The Holstein-Primakoff boson mapping |34| corresponds to the following canonical com- 
plex coordinates: 

c* = V2nv-, c = V2nv\ (4.1) 

Indeed in these coordinates the classical equations of motion are in the canonical form: 

{C\C} = i, {C*,H}=C*, {C,H}=C. (4.2) 
Through the quantization 

(Q^iU)^ (4-3) 

the classical SU(2) algebra, generated by the averages of the fermionic SU(2) algebra on the 
coherent state l^), is mapped onto a boson SU(2) algebra. Multiplying the two applications 
one obtains a mapping of the fermionic algebra onto a boson SU(2) algebra. The result is 
conventionally called as the boson expansion of the initial fermionic generator operators. 
Since this boson expansion for the generators of the SU(2) algebra has been found first by 
Holstein and Primakoff we shall refer to it as to the HP boson expansion. By the mapping 
specified by eq.(4.3), any function of the complex coordinates C*, C can be transformed into 
a function of the bosons B\ B. In particular the classical energy 7i, the average of H (2.1) 
on 1^2)5 is transformed in the HP boson expansion of the model Hamiltonian: 

H% P) = 2eB^B + Xt(B^B - ^B +2 B 2 ) 
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In the limit of f2 going to infinity the Hamiltonian H^ 1 ^ goes to: 



H ( B HP) = (2e + X 1 )B^B + \ 2 {B^ + B 2 ). 



(4.5) 



We shall refer to H^ 1 ^ as the Hamiltonian in the zeroth order boson expansion. The 
Hamiltonian in the first order boson expansion is: 

H<& = (2c + XJ&B + A 2 ^(fit 2 + 5 2 ) - ±(B* 2 N + iV5 2 ) - £-*t V. (4.6) 



4fi 



2fi 



2fi 



The transition operator linking the states described by the boson Hamiltonians (4.4), (4.5) 
and (4.6) has the expressions 



2VL 
2SL 



BW . B^B 

U p V n \l-—B + V p U n B\l 



2tt 



2VL 



U p V n B + V p U n B* 

UpVn{1 ~ ~w )B + w 5t a - 



(4.7) 

(4.8) 
(4.9) 



respectively. The matrix elements involved in the transition amplitude can be easily cal- 
culated. To save space here we describe only the results for the first order expansion. 
Diagonalizing the matrix associated to the boson Hamiltonian H BP , given in Appendix A, 
in the basis 



\m) 



'mi 



(4.10) 



for the case of mother (m) and daughter nuclei (d), one obtains the eigenstates: 



K) P = Y, C $\l), P = m,d, k = 0,1, 2, .. 



(4.11) 



Here |0) 6 denotes the vacuum state for the boson B. The corresponding eigenvalues will be 
denoted by E^f 1 . The label k = corresponds to the ground state and to simplify notations 
we shall omit it. To be more suggestive sometimes the label "0" is replaced by "g". The 
overlap of the states describing the mother and daughter nuclei have the form: 



(4.12) 
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The matrix elements of the transition operator can be easily calculated once one knows how 
to calculate the matrix elements corresponding to the basis states. These are given explicitly 
in Appendix A. The energy denominator in Eq. (2.2), corresponding to the boson expansion 
treatment, is defined as follows: 

E k = E { k m) -E^\ (4.13) 

The energy shift AE was taken equal to lMeV for small values of the proton-neutron pairing 
strength k' while for large k', half of the difference between the static ground state energies 
of mother and daughter nuclei is added to this value. The double beta transition amplitude, 
calculated by using the eigenstates of H^ P ^ is shown in the lower plot of Fig.3, lowest 
panel. Since the ground states of the mother and daughter nuclei contain many boson 
components the transition operator, although it is linear in bosons, may excite the ground 
states to a many boson state describing the intermediate odd-odd nucleus. Therefore for the 
intermediate states in Eq. (3.2) the complete set produced by the diagonalization procedure 
was considered. We found that the first excited state of the above mentioned complete set 
contributes the largest part to the double beta transition amplitude. The general trends 
shown in the upper and lower panel are similar to each other. While the left branches 
are identical, the second branch corresponding to the zeroth boson expanded Hamiltonian 
is much larger than that produced by the harmonic approximation defined by using the 
coherent state for the SU (2) group, as variational state. This indicates that for these values 
of k' the anharmonic effects prevail over the harmonic ones. While in the upper panel we 
notice an interval where the transition amplitude is not defined, the transition amplitude 
in the lower panel is vanishing in this interval. The reason is that the mixture of the 
higher boson states is large and different components have different phases and therefore 
their partial contributions cancel each other. The results for the double beta transition 
amplitude corresponding to the full HP boson expanded Hamiltonian are compared in Fig. 
4 with those obtained with the first order boson expanded Hamiltonian . One remarks that 
sizable differences are noticed only for large values of k'. It is worth to mention that while 
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in the lower panel of Fig. 3 the window where the transition amplitude vanishes is for a 



finite range of k', here this is reduced just to a point. As we showed in ref. |26|] the HP 
boson expansion is defined in that interval of k' where the standard pnQRPA does not work. 
However the final boson Hamiltonian can be diagonalized also for k' belonging to the region 
where the standard pnQRPA properly works. In this region however the boson cannot be 
interpreted as describing small oscillations around a "deformed" minimum but around the 
origin of the phase space. We would like to mention the fact that the results obtained with 
the eigenstates of the full HP Hamiltonian coincide with the exact result corresponding to 
the exact eigenstates of the fermionic quasiparticle Hamiltonian |26| . Comparing the results 
shown in Figs. 3 (first panel) and 4 one may say that in the first interval of k' the pnQRPA 
works quite well, the corresponding results lie close to the exact ones. By contrast for large 
values of k', the results corresponding to the harmonic approximation are quite different from 
the exact ones. The exact result indicate that a big effect is caused by anharmonicities. 

B. Schwinger boson expansion 



In ref. |Tj]we proposed a boson expansion formalism for treating the Gamow- Teller double 
beta transition. The underlying idea was that the boson series associated to two quasiparticle 
proton neutron operators should be chosen in such a way that the mutual commutation 
relations for the fermionic operators are satisfied in each order of the approximation. This 
condition could not be satisfied if we restricted the space to the pnQRPA bosons. However, 
the condition is fulfilled if the space is extended by including also the quadrupole charge 
conserving QRPA bosons. Therefore the boson representation of the bifermionic proton- 
neutron operators contains two types of bosons and therefore we call it Schwinger like boson 



expansion. Schwinger |35| was the first who achieved a boson mapping of the SU(2) algebra 
by using two bosons. A weak point of boson expansions for many body operators is that 
there is no quantitative measure for the "rest" when the infinite series is truncated. However, 
in the present case the evaluation of the "rest" is possible since one knows the exact result. 
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Therefore it is worth to formulate the boson expansion, used in ref. |L4| for GT double 
beta transition, for the case of Fermi transition. In this case the charge conserving boson, 
which should be added to the proton-neutron QRPA boson in order to fulfill the condition 
of preserving the mutual commutation relations, is the pairing vibrational mode. This mode 
is a particle-particle like excitation and its amplitudes are related to the overlap of the 
function describing the given nucleus with a one associated to a system obtained by adding 
(or subtracting) a pair of protons or a pair of neutrons. It is worth to remark that the 
daughter nucleus involved in a Fermi double beta decay is a component of a double phonon 
vibrational state obtained by removing a pair of neutrons and adding a pair of protons. 
Since the daughter nucleus could be reached by transforming the mother nucleus in two 
different ways, the double beta decay and the double phonon excitation, one may expect 
that the anharmonic part of the boson expanded transition operator for single beta decay 
leg might excite two boson and three boson states in the intermediate odd-odd nucleus 
having one and two charge conserving phonon factors, respectively. The pairing vibrational 
modes are usually defined by the proton-proton and neutron-neutron pairing interactions. 
Therefore there are proton and neutron vibrational modes which are decoupled from each 
other. Within the QRPA approach there are two spurious states due to the conservation of 
both the number of protons and the number of neutrons. When the space of single particle 
states is restricted to a single shell there is no physical solution. In order to allow a non- 
vanishing energy, one needs at least an additional level where to promote a pair of particles. 
To keep simplicity we maintain the single j picture but consider the contribution of the 
proton-neutron particle-hole and particle-particle interactions to the equations of motion 
of the charge conserving operators. In this way, for example, the single j of the neutron 
system will play the role of the second level for the equations of motion of the proton-proton 
monopole operators. Therefore we keep, in the quasiparticle representations of the model 
Hamiltonian, the terms with four quasiparticle operators which contribute to the equations 
of motion of two quasiparticle monopole operators: 
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H„ 



E' N + E' N 



G r 



2Vt 



(K + Vi)A\ v A w - U 2 p V p \A\ p + Kp) + ^X B 1 B pp_ 
(U A n + V*)Al n A nn - U 2 n V*(At n + Al n ) + AUXBl n B,, 
+X 1 Al n A pn + A 2 [AH + A 2 pn - Bf n - B 2 pn ] , 



-^2fi 
4 



(4.14) 



where the following notations have been used: 

Gr, 



E' p = E p + -fV p A + 2{U 2 p - V 2 )( X U 2 - XiV 2 ) 
E' n = E n + " W 2 {U 2 n - V T 2 )(x + 



E p 
A p 

V„ 

aL 



G p Vt _ G n Q 



G„nU p Vp, A n — G n QU n V r 



'Nn 

2Vt 



1 



, v p 



- alm a ~ ) T 

■ ± — t Tin rm i 

J m 



n ^->n^ L ^n ' ni 

N~ 



p,n. 




(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 

(4.20) 



Assuming quasiboson commutation relations for the operators A\ T ,A\ T one obtains the 
following linearized equations of motion: 



H 4t 

11 qj ■ ri - pp 



^(e' p -e p ) + §)aI p -( x + Xi )^aI 



T-i 7-1 " ~nn 

& p tj n 



+-^A PP -( X + Xi) 



A n A 



p^n 



E p E n 



(4.21) 



H 4t 

11 q , s\ nn 



A„A r 



ix + Xi)^/4p + I 2 « -EJ + irl 4, 



J p ± ^n 



E„ 



[H q , A 



i>p\ 



A A A 2 

^p-^n i-'n 

-fAl + ( x + Xl )^Al n 

( A 2 \ A A 

- 2{E' -13,) + -*) A pp + ( X + 



(4.22) 



Bp, 

[H q ,A nn ]=( x + Xl )^Al p -^Al n 

+(x + Xi)^Ap P -(2(E' n -E n ) + 



E p E n 



(4.23) 



En 



A r . 



(4.24) 



These equations allows us to determine the operator 
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pt - Y 4t x y 4f _ V 4 _ V 4 

tt ^-p^-pp T ^"^nn 1 p- n -pp 1 n^-nn- 



(4.25) 



which fulfills the restrictions 



h r f 



r r f 



(4.26) 



The first equation (4.26) provides an homogeneous system of linear equations for the ampli- 
tudes X and Y. To be able to solve these equations, the determinant of the coefficients has 
to be zero. This determines the excitation energy uj: 



4(4 - E p f - u? + 4^(E> - E p ) 



A(E' n -E n ) 2 -uj 2 + ^{E' n -E,, 



E„ 



-16(4 " E P )( E 'n - E n)(x + Xl) 2 ^ 

p n 



0. 



(4.27) 



The amplitudes are obtained up to a multiplicative constant which is fixed by the normal- 
ization restriction given by the second eq. (4.26 ), which reads: 



^xi + xi-y;-y^) = i. 



(4.28) 



In a similar way one determines the QRPA equations for the amplitudes and energy of the 
proton-neutron phonon operator: 



pt - v 4t _ v A 

L pn yl pn^- pn 1 pn^-pn- 



(4.29) 



The energy has the expression 



UJ 



(E' p + E' n + Ax) 2 - 4A 2 



(4.30) 



The corresponding phonon amplitudes are: 



X 



pn 



E'p + E^ + Xr- uj 



Y, 



pn 



4AI 

E' p + E' n + Ai - uj 
2A^ 



(4.31) 



Following the procedure described in detail in ref. [Ill] one obtains the following boson 
expansion of Schwinger type for the operators involved in the transition operator: 



16 



Kn = br U r l + crj„r rr + drlr pn + er pn r TT , 

A \ n = «iorJ„ + a ir pn + a 30 rj n r^r+ r + a 2 irj„r^r rr + a 12 rj„r rr r rr 

+ a 2 -iit r r^ r r pn + ai- 2 r^ T r p „r rr + ao 3 rp n r rr r rr . (4.32) 

The expansion coefficients are given in Appendix B. These boson expansion induces anhar- 
monic components for the transition operator which excite the ground state f the mother 
nucleus to some many boson states. Note that within the standard pnQRPA approach, only 
the one phonon state is accepted as intermediate state, all other states produce vanishing 
matrix elements. We denote by Mp \ k=l,2,3 the transition amplitudes for the double beta 
Fermi transition determined by one, two and three boson states of the intermediate nucleus. 
Within the first order boson expansion the result for the total amplitude is: 

M F = + M { p ] + M ( p ] , (4.33) 

where the partial amplitudes are given explicitly in Appendix C. 

Within the standard pnQRPA formalism, the double beta transitions to excited states 
of the daughter nuclei are forbidden. If the first order boson expansion of the transition 
operator is considered then the transition to the first and double phonon states are allowed. 
Here we study the transition to the double phonon pairing vibration. The amplitude for 
this transition has the expression: 

M 2 Ph = + u^OQ^Aujvy^ + u d X4i) 

I 2-/2^ ( ^ m(C3 + U ™ V n m <0) O ^ U n V p< + Wl) (4 34) 



uj™ n + 2u m + AE) 



The overlap matrices O^ 1 ^ and are defined in Appendix C. The superscripts m and d 
suggest that the given quantities characterize the mother and daughter nuclei, respectively. 

The results corresponding to transition operators defined in the standard pnQRPA (i.e. 
linear in the pn bosons), and higher order pnQRPA by including quadratic and cubic terms 
in the bosons expansion are shown in Fig. 5 together with the exact result. Of course since 
the pn boson collapses at about k' = l.OBMeV this boson expansion is justified only for 
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the interval below this value. One notes that the quadratic terms in bosons included in the 
transition operator modifies the pnQRPA result for the double beta transition amplitude 
such that the final result is in excellent agreement with the exact one. Also one notices that 
the third order boson terms are practically not modifying the result obtained with the linear 
and quadratic terms. One may conclude that at least for the single j shell considered here 
and the double beta Fermi transition the first order boson expansion is approaching very 
well the exact result. Of course important deviations appear when we approach the critical 
value of k' where the pn boson collapses. For example for k' where the pnQRPA amplitude is 
equal to zero the exact result is about 0.9 MeV -1 . As we have already mentioned within the 
pnQRPA approach the transition to an excited state of the daughter nucleus is forbidden. 
However this is allowed within the boson expansion formalism. Here we calculated the 
amplitude for the transition Of — > (2ph) f from the ground state of the mother to the two 
vibrational phonon states of the daughter. Such a state may be alternatively reached in a 
process which removes a pair of neutrons and adds a pair of protons to the mother nucleus. 
Therefore if pairing vibrational phonons play a role by inducing anharmonic effects in the (3~ 
transition operator the double phonon state might be populated in a double beta decay. The 
results are shown in Fig 6. From there one remarks that the process amplitude is reasonable 
large and if that order of magnitude persists for realistic calculations one might hope that 
such transitions can be experimentally identified. Coming back to Fig 5, one notices that 
around k' — 1.2 the exact amplitude suggests a phase transition since its first derivative has 
a discontinuity there. In the RPA procedure such a transition is associated to a Goldstone 
mode which appears whenever a new symmetry open up. What happens however in the 
case of the exact solution? To give an answer to this question we plotted in Fig 7 the first 
five eigenvalues of the HP boson expanded Hamiltonian for the mother (upper panel) and 
the daughter (lower panel) nuclei. Also in Fig. 8 we give the first two leading amplitudes 
of the first four eigenstates. In Fig. 7 one sees that not far from the turning point in Fig. 
5 the first excited state is almost degenerate with the ground state. The virtual transition 
to such a state is almost forbidden, despite the fact that its excitation energy is so small, 
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since the matrix elements of the transition operator are suppressed. This happens due to 
the fact that two amplitudes of the wave functions (see Fig. 8) in this region are becoming 
equal in magnitude but of opposite phase. Indeed for a small k' the dominant component 
of the ground state is the vacuum state while for k' greater than 1.75 the third component, 
corresponding to the two phonon states, becomes dominant. That is in fact a signature for 
a super-deformed ground state. A similar picture holds also for the first excited state. For 
small k' this is a one boson state while for k' larger than 1.5 it is a three boson state. It is 
interesting to see that for the mother nucleus the initial components ordering of the third 
excited state is recovered by increasing k' beyond the value 2.4. It seems that beyond the 
k'=1.2 the contributions of the higher boson states to the total transition amplitude are 
constructively interfering. 

Let us now conclude the present Section. While in the preceding section we treated the 
double beta transition amplitude within an QRPA like approaches provided by VP1 and 
VP2, here boson expansion procedures have been used in order to incorporate higher RPA 
correlations. In the region of large particle-particle interaction strength k' we have used a full 
HP boson expanded Hamiltonian and two truncated HP boson expansions. In the interval of 
small strengths for the particle-particle interaction k' a Schwinger like boson expansion have 
been used. In both cases we compared the results of truncated boson expansions with the 
corresponding RPA like approach as well as with the exact results. In both cases the first 
order boson expansions for the transition operator provides a description which is very close 
to the exact picture. While in the interval of small particle-particle interaction strength k' 
the VP1 result is not bad comparing it to the exact one, in the complementary region the 
result produced by the VP2 exhibits large deviation from the exact result. This reflects the 
fact that the exact result includes large contributions coming from anharmonicities. The 
transition from one regime to another could be interpreted by inspecting the structure of 
the first exact eigenstates of the model Hamiltonian. Indeed going beyond k' 1.5 the 
dominant component in the ground state is the two phonon state while in the first excited 
state the three phonon state prevails. 
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V. NEW FEATURES OF THE RPA VACUA AND STATIC GROUND STATES 



The BCS approximation determines variationally the ground state of a system of nu- 
cleons interacting among themselves through pairing force. This approximation determines 
simultaneously the optimal ground state and a new type of particle excitations which admit 
the found ground state as vacuum . In a time dependent variational formalism the BCS 
ground state turns out to be the static ground state with respect to which one could de- 
fine the small oscillations which account for the quasiparticle two body interaction. This 
approach is fully equivalent to the QRPA formalism. The QRPA formalism finds also an 
excitation operator which involves excitations of many quasiparticles from the ground state. 
However the vacuum state, or in other words the new ground state involving correlations 
corresponding to these degrees of freedom, is not a solution of a variational equation. 

In ref. ||26|| , using 4 different trial functions we determined variationally four static ground 



states and moreover four RPA like phonon states. The vacuum states for the phonon opera- 
tors yielded by VP1 and VP2 can be analytically obtained. Indeed the equation expressing 
the condition of being vacuum for the phonon operators can be solved and the results are: 

Y 



\0) w = A/Wexp 

|0)au(2) = A/" sn ( 2 )exp 



w f /it \2 

2X„ l pn) . 



2X, - pn 



(5.1) 
(5.2) 



L s«(2) 

The indices w and su(2) suggest that the amplitudes X and Y characterize the phonon 
operator describing small oscillations around a static ground state represented by a coherent 
state of the Weyl and SU(2) groups, respectively. These amplitudes were determined in our 
previous publication (see eq. 6.20 and 6.28 of ref. |26|| ). Assuming for the operators 
involved in (5.1,5.2) quasiboson commutation relations one obtains the following expressions 
for the norms. 

n y on 

^ V 9/1 

^( 2 )- 2 = E( 2 ^) 2/ 7=^ 5 (5-3) 

; =0 2A su(2 ) [ll) 
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It is worth noting that both vacua are exponential functions of (A^) 2 . On the other hand 
so are the static ground states yielded by the VP3 and VP4 approaches: 

o 

|*s) = Ku(i,i) exp a (Ajj 2 , a = —tanh(2p) (5.4) 
|tf 4 > =A/" 4 expz(4J 2 , z = pe i +. (5.5) 

The norm A/" s «(i,i) has been analytically calculated 

A/; u( i,i) = , , (5.6) 

Jcosh(2p) 

while A/4 has only been numerically determined. The parameters entering the defining 
equations (5.4) are those which produce a minimum energy for our system. The question 
we address in this paper is how one compares the vacua of the phonons defined by means 
of the VP1 and VP2 formalisms to the static ground states of VP3 and VP4. If they are 
close to each other one could state that the solutions of the variational equations VP3 and 
VP4, i.e. the static ground states, approximate the vacua of VP1 and VP2 which are the 
dynamic ground state of the system. For an easier presentation it is convenient to use the 
unified notation 

|$ fe ) =Af k exp[a k (Alf] ,fc = l,..,4. (5.7) 

The indices k — 1,2 correspond to the vacua defined within the VP1 and VP2 formalisms, 
respectively, while k — 3,4 give the static ground states provided by VP3 and VP4, respec- 
tively. In Figs 9,10 and 11 we compare the norms, exponents and the products Afa of the 
four functions. Actually the first and last quantities are just the weights of quasiparticle 
vacuum and double phonon excitations, respectively, for the four states mentioned above. 
From these figures one notes that for situations when the first two components of these 
four functions are dominant one may state the following. For k' < 0.8MeV the vacuum 
state defined by the VP1 approach is quite well approximated by the static ground state 
corresponding to VP3. The static ground state of the VP4 approach approximates very well 
the vacuum state of VP1 and reasonable well the one corresponding to VP2 in the interval 
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of k' from 1.3 to 1.8 MeV. This result is very important since one could think of formulating 
a theory which improves the RPA approach by that the vacuum of the phonon operator 
satisfy the equations provided by a variational principle. 

VI. CONCLUSIONS 

The main results obtained in the previous sections can be summarized as follows. The 
model Hamiltonian for a single j shell was treated within a harmonic picture for values of 
the particle-particle interaction strength k' below and beyond the critical value, where the 
standard pnQRPA ceases to be any longer valid. The energies and the states describing the 
proton and neutron systems are obtained through a time dependent variational principle 
choosing as trial function coherent states of the symmetry groups of Weyl and S77(2), re- 
spectively. The results are used to calculate the transition amplitude for double beta Fermi 
transitions. In the range for larger k' values one may define a Holstein-Primakoff boson 
representation of the model Hamiltonian which is further on abusively used also for the 
complementary interval. The transition amplitude for the double beta decay was calculated 
using the zeroth, the first order as well as the full HP boson expansion for the model Hamil- 
tonian. The last procedure reproduces the exact result. For the first part of interval we also 
considered a Schwinger type boson expansion which is similar to the one proposed for the 
treatment of the Gamow Teller double beta decay, by two of the present authors (A. A. R. 
and A. F.). By the boson expansion analysis we concluded that the first order truncation 
provides a good approximation to the exact results. In the case of the Schwinger boson 
expansion we calculated the transition to the double phonon vibrational state. Since the 
daughter nucleus can be reached by a double phonon excitation of the mother nucleus, the 
vibrational states should play an important role in explaining quantitatively the decay rates. 
The two harmonic approximations, for small and large values of the particle-particle two 
body interaction strength k', define phonon operators whose vacua are dynamical ground 
states. These states can be analytically expressed and they exhibit a striking resemblance 
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with the variational states used by the VP3 and VP4 approaches. The question is how one 
compares these states with each other. We identified intervals where the dynamic ground 
states are approximated quite well by the states determined variationally. This is in fact 
the first step toward solving an old standing problem of determining variationally the RPA 
ground state. We hope that the results of this paper will stimulate further studies both for 
the field of double beta decay in the region of large particle-particle interaction but also for 
the static and dynamic ground states of many body systems. 

VII. APPENDIX A 



The matrix elements of the first order expanded Hamiltonian in the many boson basis 
have the following explicit expressions: 



(m\Hffp\m) = [2e + \\)m — ^m(m — 1), 



Ai 



(m\Hjp P \m + 2) = A 2 4 ^ 1 y/ (m + l)(m + 2) - ^mym + l)(m + 2), 



(m + 2\H { H ) P \m) = (m\H ( H P \m + 2) 



(i) 



(7.1) 



The matrix elements for the first leg and second leg of the double beta transition are given 
by: 



1,1' 

1TYY J71 1 

u P v n (i - —)VTTT6,, l+1 + v p u n {\ - - 1?r 



)\Tl5 V) i-i 



d <oi|/?t|ot> d = v^£0? ( 



(d) n (d) 
ki 



i,i' 



x 



u p v n (i - ^)VTTTs l/;l+1 + v p u n (i - ^1)^/16^ 



(7.2) 



VIII. APPENDIX B 

Here we list the coefficients defining the boson expansions (4.32). 
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A . (^pYpn 4~ Y fl X ptl) . 

J 

2 

~(X n Xp n -\- YpYpn) ; 

2 

~ (XpXp n -\- YnYpn)i 
J 

2 

~(X n Yp n "I - ^p^pn) ; 



^10 — -X"; 



pn; 



&01 — 



2X. 



a-30 

«21 
Ol2 
«21 
«12 
003 



prt 



2X 



p« 



2X 



ph 



2K 



(X p Yp + X n Y^), 
(Xp + X 2 n + Y p 2 + Y% 
(XpYp + X n Y n ), 



n 

2K 



(XpYp -\- X n Y n 



n 



' m (Xp + Xl + Yp + F„ 2 ) 



2Y 



n 



(XpYp -\- X n Y n 



(8.1) 



In what follows we shall briefly describe the procedure presented in ref. [pL4|,p~5|| . From the 
eqs.. (4.32) one derives the equations: 



b 

a 30 



<0| 
1 



b 1 r 



10), 



(0| 



r 4t 

- 1 - p?l; ^pn 



|0>. 



(8.2) 



In the above equations, all commutators are exactly evaluated except for the last one for 
which the quasi-boson approximation is used. It is worth noting that the coefficients cal- 
culated in this way determine a boson representation for the bi-fermionic operators which 
satisfy the mutual commutation relations in the first order. 



IX. APPENDIX C 



Here we list the transition amplitudes corresponding to the pnQRPA and first order 
boson expansion: 
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(TT m \f m Y m _|_ rymT/m /jjd\rdvd i Tjd^id^d \ 

1 r,,^ \^ n " p pn ' ^ p *n pn)^md\^n * p pn ' ^ p " n pn) 



Ml = 2tt 



Mi 2) = 2VL- 



u m + AE 

(2), 



_ V rn V rn e m + U^U^^O^-V^V 1 + U d p U^) 



u™+u m + AE 



M ( 3) = lcl (u?v?<® + U£}T<}o^(uX< + ujvM fg n 



u;™ + 2cu m + AE 



The overlap matrices have the expressions: 



(1^ = ( X m X d — Y m Y d \ 

^md \ pn pn pn pn/i 

0(2) ym | \rm \rd \rm\fd 

md - A p A p + A n A n _ Y p Y p ~ Y n Y n > 

0(3) / \rm vd -\/~mxrd \ ( vm vd . vmyrf v^ra-iAd v m vd\2 m o\ 

- l A pn A pn _ r P n VJl A p A p + A n A n ~ Y p Y p ~ Y n Y n ) ■ 
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FIGURES 




FIG. 1. The energies of the harmonic modes corresponding to the VP1 (left branch) and VP2 
(right branch) defined in eqs. (2.7) and (2.8) are plotted as functions of the proton-neutron pairing 
strength (eqs. (2.2) and (2.11)) k' for the mother (a) and the daughter (b) nuclei. Energies as well 
as strength x' are given in units of MeV 
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FIGURE 2 



k=3 according to eq. (5.4) (x'=0.0) 




FIG. 2. Energies of the harmonic modes corresponding to VP3 (point and dashed-point lines) 
and VP4 (full and dashed lines) are given in units of MeV as functions of proton-neutron pairing 
strength k' 
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FIGURE 3 
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FIG. 3. The double beta Fermi transition amplitudes obtained with the initial, intermediate 
and final states provided by VP1 (full line) VP2 (dashed line) formalism (a) and by diagonalizing 
the limit Q goes to infinity of the HP boson expanded Hamiltonian (b) are given as function of the 
particle-particle interaction strength k'. 
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FIGURE 4 




0.5 1 1.5 2 

k' [MeV] 

FIG. 4. The double beta transition amplitude obtained by considering the states for mother, 
daughter and intermediate odd-odd nuclei as eigenstates of the Holstein-Primakoff boson expanded 
Hamiltonian (full line) and of the first order truncated HP Hamiltonian (dashed line) as function 
of the proton-neutron pairing strength k' defined as in eq. (2.11) 
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FIGURE 5 
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FIG. 5. The transition amplitude for double beta Fermi decay calculated with the standard 
pnQRPA (long dashed line), pnQRPA plus the second order corrections (dashed line) and pnQRPA 
plus second and third order corrections (dashed-point-point line)in the transition operator is rep- 
resented as function of k'. For comparison the exact result is also presented (full line). 
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FIGURE 6 
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FIG. 6. The transition amplitude from the ground state of the mother nucleus to the two 
pairing vibrational states, characterizing the daughter nucleus, is plotted as a function of the 
proton- neutron pairing strength k'. 
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FIGURE 7 




FIG. 7. The first five eigenvalues of the Holstein-Primakoff boson expanded Hamiltonian (4.4), 
for mother (a) and daughter (b) nuclei, are represented as functions of k'. 
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FIGURE 8 
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FIG. 8. The absolute values for the leading two amplitudes of the first four eigenstates of the 
Holstein-Primakoff boson expanded Hamiltonian given by Eq. (4.4) are represented as functions 
of k'. In the upper panel the 1st and 3rd eigenstates are analyzed while in the lower panel the 2nd 
and the 4th are given. 
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FIGURE 9 
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FIG. 9. The norms of the states defined by eqs.. (5.1), (5. 2), (5. 4), (5. 5) are plotted as functions 
of the proton- neutron pairing strength k'. 



36 



FIGURE 10 
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FIG. 10. The exponent cc*. entering the definition of the states \$k) given by Eq. (5.7). 
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FIGURE 11 
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FIG. 11. The weights for a two phonon state in the structure of the states are plotted as 
function of the strength of the particle particle interaction. 
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